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INTRODUCTION 

This chapter is devoted to exact conditions in time- 
dependent density functional theory. Many conditions 
have been derived for the exact ground-state density 
functional, and several have played crucial roles in the 
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construction of popular approximations. We believe that 
the reliability of the most fundamental approximation of 
any density functional theory, the local density approxi- 
mation (LDA), is due to the exact conditions that it sat- 
isfies. Improved approximations should satisfy at least 
those conditions that LDA satisfies, plus others. (Which 
others is part of the art of functional approximation) . 

In the time-dependent case, as we shall see, the adia- 
batic LDA (ALDA) plays the same role as LDA in the 
ground-state it satisfies many exact conditions. 

But we do not have a generally applicable improvement 
beyond ALDA that includes nonlocality in time. For 
TDDFT, we have a surfeit of exact conditions, but that 
only makes finding those that are useful to impose an 
even more demanding task. 

Throughout this chapter, we give formulas for pure 
DFT for the sake of simplicity (e.g. E'xcM), but in prac- 
tice spin DFT is used (e.g. i?xc["-t, %])■ We use atomic 
units everywhere (e^ = h = rrie = 1), so energies are in 
units of Hartrees and distances are in Bohrs. 



I. REVIEW OF THE GROUND STATE 

In ground-state DFT, the unknown exchange- 
correlation energy functional, Exc[n], plays a crucial role. 
In fact, it is this energy that we typically wish to approx- 
imate with some given level of accuracy and reliability, 
and not the density itself. Using such an approximation 
in a modern Kohn-Sham ground-state DFT calculation, 
we can calculate the total energy of any configuration of 
the nuclei of the system within the Born-Oppenheimer 
approximation. In this way we can extract the bond 
lengths and angles of molecules and deduce the lowest en- 
ergy lattice structure of solids. We can also extract forces 
in simulations, and vibrational frequencies and phonons 
and bulk moduli. We can discover response properties 
to both external electric fields and magnetic fields (using 
spin DFT) . The accuracy of the self-consistent density is 
irrelevant to most of these uses. 

Given the central role of the energy, it makes sense to 
devote much effort to its study as a density functional. 
Knowledge of its behavior in various limits can be cru- 
cial to restraining and constructing accurate approxima- 
tions, and to understanding their limitations. This task 
is greatly simplified by the fact that the total ground- 
state energy satisfies the variational principle. Many ex- 
act conditions use this in their derivation. 

In this section we will review some of the more promi- 
nent exact conditions. They almost all concern the en- 
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ergy functional, which, as mentioned above, is crucial for 
good KS-DFT calculations. We also refer the interested 
reader to Ref. [1] for a thorough discussion. First, we will 
go over some of the formal definitions in DFT. 



I.l. Basic definitions 



The XC energy as a functional of the density is written 
as [2, 3] 



min (*|'f + V;c|*)-TsN-f/[n] 



(1) 



where ^' is a correctly antisymmetrized electron wave- 
function, the minimization of the kinetic and electron- 
electron repulsion energies is done over all such wavefunc- 
tions that yield the density n{r), Ts[n] is the minimum 
(non-interacting) kinetic energy of a system with density 
n(r), and 



U[n] 



3 , n[r)n{v') 

d r — — 



(2) 



is the Hartree energy. The XC energy is usually split 
into an exchange piece, i?x, and a correlation piece, = 
Exc — E^. Exchange can be defined in a HF-like way in 
terms of the KS spin orbitals (/)icr(r): 



Ex 



Ir-r'l 



(3) 

To perform the self-consistent calculations in the non- 
interacting system, we need the functional derivative of 
the XC energy. 



WxcW(r) = 



5Exo[n] 
5n{r) 



(4) 



This is called the XC potential, and it is the essential 
part of the multiplicative KS potential fs['T.](r). 

Orbital dependent functionals: Some functionals 
are most naturally expressed in terms of the orbitals 
rather than the density. When varying the orbitals of 
these functionals, nonlocal potentials are obtained. For 
example, varying (pi^ in Eq. (3) leads to the nonlocal ex- 
change term used in HF. There is a way to transform such 
orbital-dependent functionals into local potentials as in 
Eq. (4). This procedure is known as optimized effective 
potential (OFF) or optimized potential method (0PM) 
and is computationally expensive [4]. Using OFF for E^ 
results in the exact exchange approximation (EXX) for 
Exc in KS-DFT. The Krieger, Li, and lafrate (KLI) ap- 
proximation is a way to approximately solve EXX [5]. 

Adiabatic connection: One can imagine smoothly 
connecting the interacting and non-interacting systems 
by multiplying the electron-electron repulsion term by 
A, called the coupling-constant. Changing A varies the 



strength of the interaction, and if we simultaneously 
change the external potential to keep the density fixed, 
we have a family of solutions for various interaction 
strengths. This makes all quantities (besides the density) 
functions of A. When A = 0, one has the non-interacting 
KS system, and when A = 1, one has the fully interacting 
system. The following coupling-constant relations hold. 

XC energy A dependence: Altering the coupling- 
constant is simply related to scaling the density: 

E^M^^^Exa[ni/x], (5) 
where ni/xir) is the scaled density 

"7(1") = 7^ "(71") > (6) 

with 7 = 1/A. 

Adiabatic connection formula: By using the 
Hellmann-Feynman theorem, one can show: 



Exo[n] = / dAC/^Jn]/A 



(7) 



where U^^ is the potential contribution to exchange- 
correlation energy (J7xc = Vcc — U) at coupling-constant 
A. 



1.2. Standard approximations 

Despite a plethora of approximations [6], no present- 
day approximation satisfies all the conditions mentioned 
in this chapter, as seen in tests on bulk solids and surfaces 
[7]. With that the case, one must choose which conditions 
to impose on a given approximate form. Non-empirical 
(ab initio) approaches attempt to fix all parameters via 
exact conditions [8, 9], while good empirical approaches 
might include one or two parameters that are fit to some 
data set [10-12]. 

There are two basic flavors of approximations: pure 
density functionals, which are often designed to meet con- 
ditions on the uniform gas, and orbital-dependent func- 
tionals [13], which meet the finite-system conditions more 
naturally. The most sophisticated approximations being 
developed today use both [14]. For a good discussion on 
what approximation is the right tool for the job, see Ref. 
[15]. 

LDA: The local density approximation is the bread 
and butter of DFT. It is the simplest, being derived from 
conditions on the uniform gas [16]. Though it is too 
inaccurate for quantum chemistry (being off by about 
1 eV or 30 kcal/mol), it is useful in solids and other bulk 
materials where the electrons almost look like a uniform 
gas. There can only be one LDA. 

GGA: The generalized gradient approximation came 
from trial and error when energies were allowed to depend 
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on the gradient of the density. While more accurate than 
the LDA (getting errors down to 5 or 6 kcal/mol), and 
thus useful for quantum chemistry applications, there is 
no uniquely-defined GGA. BLYP is an empirical GGA 
that was designed to minimize the error in a particular 
data set. PBE is a non-empirical GGA designed to satisfy 
exact conditions. 

Hybrid: Hybrids have an exchange energy which is a 
mixture of GGA and HF, which attempts to get the best 
of both worlds: 

E^l^^E^^^ + aiE,^E^^% (8) 

where E^ is defined in (3). The parameter a was argued 
to be 0.25 for the non-empirical PBEO, but is fitted for 
the empirical B3LYP. 



1.3. Finite systems 

The following conditions are derived for finite systems, 
just as the Hohenberg-Kohn theorem is. 

Signs of energy components: From the variational 
principle and other elementary considerations, one can 
deduce 

Ey,c[n] < 0, Ea[n] < 0, < 0. (9) 

Zero XC force and torque theorem: The XC po- 
tential cannot exert a net force or torque on the electrons 
[17]: 

j d^r n{r) Vi;xc(r) = 
J d^r n(r) r x Vwxc(r) = 0. (10) 



is in the high density limit (7 large) or low density limit 
(7 small) [17, 18]: 

^c["7] <7^c["] (7<1) 
Eo[n^]>lEo[n] (7 > 1) 

E^ [n^] = [n] + 4'^ ["]/7 + • • • (7 ^ ^) 
Ec[n^]^^B[n]+^^/^C[n] + --- (7^0), (13) 

where Eiy[n], El?'[n], B[n], and C [n] are all scale- 
invariant functionals. These conditions are depicted in 
Fig. 1. Not all popular approximations satisfy these con- 
ditions. 

Correlation Energy Scaling 
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FIG. 1. Scaling of the correlation energy in ground state 
DFT, as well as the various conditions from Eq. (13). The 
first two relations are illustrated with the dotted line. For 
7 < 1, the exact curve (solid) must lie below this dotted line, 
and for 7 > 1 the exact curve must lie above - in both cases 
within the shaded region of the graph. The high density limit 
is shown with the dot-dashed line, and the low density limit 
with the dashed line. It is believed that not only is Ec[n-y] 
monotonic, but also its derivative with respect to 7. Color 
online. 



XC virial theorem: 

E^c[n]+Tc[n] ^- j d'V n(r) r • Vwxc(r), (11) 

where Tc = T — Tg is the kinetic contribution to the 
correlation energy. The XC virial theorem as well as the 
zero XC force and torque theorem are satisfied by all 
sensible approximate functionals. 

Exchange scaUng: By using the scaled density (6), 
one can easily show 

E^[n^] ^ -1 E^[n]. (12) 

Correlation scaling: The scaling of correlation is less 
simple than exchange, and will depend on whether one 



Self-interaction: For any one-electron system [19], 
E^[n]^-U[n], Ec = (N = 1). (14) 

Lieb— Oxford bound: For any density [20], 

E^c[n] > 2.273 E^^^[n]. (15) 

In addition to conditions on i?xci '^^ ^Iso know some 
exact conditions on the XC potential and the KS eigen- 
values. 

Asymptotic behavior of potential: Far from a 
Coulombic system 

v^c{r) ^ -1/r (r^oo), (16) 

and 

ehomo = —I, (17) 



3 



Overview for TDDFT 



where enoMO is the position of the highest occupied KS 
molecular orbital, and / the ionization potential. These 
results are intimately related to the self-interaction of one 
electron. 



1.4. Extended systems 

The basic theorems of DFT are proven for finite quan- 
tum mechanical systems, with densities that decay at 
large distances from the center. Their extension to ex- 
tended systems, even those as simple as the uniform gas, 
requires careful thought. For ground-state properties, 
one can usually take results directly to the extended limit 
without change, but not always. For example, the high- 
density limit in Eq. (13) of the correlation energy for a 
finite system is violated by a uniform gas. With these 
things in mind, we will now discuss a set of conditions 
that involve the properties of the uniform or nearly uni- 
form electron gas. 

Uniform density: When the density is uniform, 
E^c — Cxc'^l")"'^' where e™''(") is the XC energy den- 
sity of a uniform electron gas of density n, and V is the 
volume. This forms the basis of LDA. 

Slowly varying density: For slowly varying densi- 
ties, Exc should recover the gradient expansion approxi- 
mation (GEA): 

Exc[n] = J d^reZ'Hn) + J d^r Ae'^^^{n,Vn) + ■■■ 

^El^^^[n]+AE^^^[n] + -- - , (18) 

where Ae'^^^ln, Vn) is the leading correction to the LDA 
XC energy density for a slowly varying electron gas [21]. 
However, the GEA was found to give poor results and vi- 
olate several important sum rules for the XC hole when 
applied to other systems [22] . Fixing those sum-rules led 
to the development of ab initio GGAs. Though impor- 
tant in obtaining the energy for the ground-state, the XC 
hole rules have not been used in TDDFT and therefore 
will not be further discussed in this chapter. 

Linear response of uniform gas: Another generic 
limit is when a weak perturbation is applied to a uniform 
gas, and the resulting change in energy is given by the 
static response function, xiQy'^ — 0)- This function is 
known from accurate Quantum Monte Carlo calculations 
[23], and approximations can be tested against it. 



II. OVERVIEW FOR TDDFT 

The time-dependent problem is more complex than the 
ground-state problem, making the known exact condi- 
tions more difficult to classify. We make the basic dis- 
tinction between general time-dependent perturbations. 



of arbitrary strength, and weak fields, where linear re- 
sponse applies. The former give conditions on Vxc[n]{r, t) 
for all time-dependent densities, the latter yield condi- 
tions directly on the XC kernel, which is a functional 
of the ground-state density alone. Of course, all of the 
former also yield conditions in the special case of weak 
fields. 

In the time-dependent problem, we do not have the 
energy playing a central role. Formally, the action plays 
an analogous role (see van Lccuwcn ch 6), but in prac- 
tice, we never evaluate the action in TDDFT calculations 
(and it is identically zero on the real time evolution). In 
TDDFT, our focus is truly the time-dependent density it- 
self, and so, by extension, the potential determining that 
density. Thus many of our conditions are in terms of the 
potential. 

Most pure density functionals for the ground-state 
problem produce poor approximations for the details of 
the potential. Such approximations work well only for 
quantities integrated over real space, such as the energy. 
Thus approximations that work well for ground-state en- 
ergies are sometimes very poor as adiabatic approxima- 
tions in TDDFT. Their failure to satisfy Eq. (16) leads 
to large errors in the KS energies of higher-lying orbitals 
(for example, consider the LDA potential for Helium in 
Figure 3 of Ref . [24] , which falls off exponentially rather 
than as — 1/r), and (17) is often violated by several eV. 

In place of the energy, there are a variety of physical 
properties that people wish to calculate. For example, 
quantum chemists are most often focused on the first 
few low-lying excitations, which might be crucial for de- 
termining the photochemistry of some biomolecule. Then 
the adiabatic generalization of standard ground-state ap- 
proximations is often sufficient. At the other extreme, 
people who study matter in strong laser fields are often 
focused on ionization probabilities (see Ullrich and Ban- 
drauk chapter), and there the violation of Eq. (17) makes 
explicit density approximations too crude, and requires 
orbital-dependent approximations instead. 



II. 1. Definitions 

In contrast to the ground-state problem, the XC po- 
tential depends not only on the density but on the ini- 
tial wavefunction ^'(0) and KS Slater determinant $(0), 
written symbolically as Vxc\n\ ^(0), $(0)](ri:). This more 
complicated dependence comes about because two differ- 
ent wavef unctions, which are chosen to have the same 
density for all time, can come from completely different 
external potentials, which the XC potential accounts for. 
We can get rid of this initial wavefunction dependence if 
we start from a non-degenerate ground-state, where the 
wavefunction is a functional of the density alone, via the 
Hohenberg-Kohn theorem [25] . These things are further 
discussed in Neepa's chapter. 

As the density evolves, the XC potential is determined 
not solely by the present density n(r, i), but also by the 
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history n{r,t') for < t' < t. However, it is useful to 
break the XC potential up into two pieces, an adiabatic 
piece which only deals with the present density, and a dy- 
namic piece which incorporates the memory dependence: 



v,c[n; *(o), mm - vtrin; MO), mm 



The adiabatic piece of the potential, 



,,adi; 



[n](rt) 



Sn{r) 



(19) 



(20) 



nit) 



is the XC potential for electrons as if their instantaneous 
density were a ground state. In the spirit of DFT, the 
dynamic piece is everything else. 

In the linear response regime, small enough pertur- 
bations to the density will continuously change the XC 
potential: 

Vxc[n + 5n]{rt) - v^cMirt) = 

dt' / dV/xcW(r,r';t,i')'5ri(r'i'), (21) 



where /xc is the XC kernel, which can be written formally 
as the functional derivative: 



/xcK](r,r';i,0 = 



(5wxc W(rO 



Sn{r't') 



(22) 



The evaluation at ng reminds us that /xc is used for 
the linear response of a density variation away from a 
ground-state density no. 

Like the XC potential, the kernel can also be broken 
down into an adiabatic piece: 



/^^-(r,r';t,i') = 



Sn{r) Sn{r') 



Sit~t'), (23) 



lit) 



and a dynamic piece, which includes memory and every- 
thing else. The kernel is often Fourier-transformed from 
position space in the relative coordinate (r — r') to mo- 
mentum space (with wave- vector q), from the relative 
time (t — t') to frequency (cj) domain, or both. Some 
conditions are more naturally expressed in momentum 
space and/or in the frequency domain. In the frequency 
domain, the adiabatic piece can be written as 



/x^r(r,r')= lim/xc(r,r». 



(24) 



The kernel is discussed in more detail in Chapter 4 
(TDDFT intro by Gross). 



II. 2. Approximations 



neglects the dynamic term of Eq. (19); these classes are 
respectively adiabatic and non-adiabatic (i.e. memory) 
approximations. In the adiabatic approximation, famil- 
iar ground-state functionals (such as LDA, GGA, and 
hybrids) can produce XC potentials when one uses the 
approximate i?xc in Eq. (20). We mention two notable 
adiabatic approximations now. 

ALDA: The prototype of all TDDFT approximations 
is the Adiabatic Local Density Approximation, and it is 
the simplest pure density functional. The XC potential 
is as simple as can be: 



,.ALDA 



n] (ri) 



de-"(n) 



dn 



(25) 



n(rt) 



In linear response, the ALDA kernel is 



/x^ci^°^(r,r';i,t')^ ^''"""^"^ 



d'n? 



n(rt) 



S^{r-r')5{t-t'). 

(26) 

Like its ground-state inspiration, ALDA satisfies impor- 
tant sum rules by virtue of its simplicity, namely its lo- 
cality in space and time. ALDA is commonly used in 
many calculations, and is described further in chap 1. 

AA: In the 'exact' adiabatic approximation, we use 
the exact E^c in Eq. (20). This approximation is the best 
that an adiabatic approximation can do, unless there is 
some lucky cancellation of errors. Hessler et al. [26] in- 
vestigated AA applied to a time-dependent Hooke's atom 
system and found large errors in the instantaneous cor- 
relation energy. For the double ionization of a model 
Helium atom, Thiele et al. [27] discovered that non- 
adiabatic effects were important only for high-frequency 
fields. 

A key aim of today's methodological development is 
to build in correlation memory effects. Any attempt to 
build in memory goes beyond the adiabatic approxima- 
tion, and thus belongs in the non-adiabatic class of ap- 
proximations. The next three approximations belong to 
this dynamic class. 

GK: The Gross-Kohn approximation is simply to use 
the local frequency-dependent kernel of the uniform gas. 



f^^{vy-M^5\v-v')f^^\n{v)-Lo), (27) 



where 



/xT(^;c^)^iim/xT(«;9,a.) 



(28) 



is the response of the uniform electron gas with density 
n. GK was the first approximation to go beyond the 
adiabatic approximation, but was found to violate trans- 
lational invariance. 



As we go through the various exact conditions, we will 
discuss whether the simplest approximations in present 
use satisfy them. We can divide all approximations into 
two classes based on whether or not the approximation 



VK: The Vignale-Kohn approximation sought to im- 
prove upon the shortcomings of GK. The VK approxi- 
mation is simply the gradient expansion in the current 
density for a slowly- varying gas (see Vignalc chapter). 
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XX: Exact exchange, the orbital-dependent func- 
tional, is treated as an implicit density functional (see 
Kiimmcrs orbital chapter (11)). When treated this way, 
XX has some memory for more than two unpolarized 
electrons. With the exception of XX, non-adiabatic ap- 
proximations are usually limited to the linear response 
regime and approximate the kernel, /xc- There is now 
a major push to go beyond linear response for non- 
adiabatic approximations. The first such attempt was 
a bootstrap approach of Ref. [28] . More recent attempts 
are described in Chapter 26 (Tokatly) of the book and in 
Rcf. [29]. 

III. GENERAL CONDITIONS 

In this section, we discuss conditions that apply no 
matter how strong or how weak the time-dependent po- 
tential is. They apply to anything: weak fields, strong 
laser pulses, and everything in between. They apply also 
to the linear response regime, yielding the more specific 
conditions discussed in Section IV. 



III.l. Adiabatic limit 

One of the simplest exact conditions in TDDFT is the 
adiabatic limit. For any finite system, or an extended 
system with a finite gap, the deviation from the instan- 
taneous ground-state during a perturbation (of arbitrary 
strength) can be made arbitrarily small. This is the 
adiabatic theorem of quantum mechanics, which can be 
proven by slowing down the time-evolution, i.e., if the 
perturbation is V{t), replacing it by V{t/T) and making 
T sufficiently large. 

Similarly, as the time-dependence becomes very slow 
(or equivalently, as the frequency becomes small), for 
such systems the functionals reduce to their ground-state 
counterparts: 

vxc (r , -> v^c Ht)] (r) (r ^ oo) (29) 

where ■i;xc["'](r) is the exact ground-state XC potential 
of density n{r). 

By definition, any adiabatic approximation satisfies 
this theorem, and so does XX, by reducing to its ground- 
state analog for slow variations. On the other hand, if 
an approximation to Vxci^t) were devised that was not 
based on ground-state DPT, this theorem can be used 
in reverse to define the corresponding ground-state func- 
tional. 



III. 2. Equations of motion 

In this section, we discuss some elementary conditions 
that any reasonable TDDFT approximation should sat- 
isfy. Because these conditions are satisfied by almost all 



approximations, they are best applied to test the qual- 
ity of propagation schemes. For a scheme that does not 
automatically satisfy a given condition, then a numeri- 
cal check of its error provides a test of the accuracy of 
the solution. A simple analog is the check of the virial 
theorem in ground-state DFT in a finite basis. 

These conditions are all found via a very simple proce- 
dure. They begin with some operator that depends only 
on the time-dependent density, such as the total force on 
the electrons. The equation of motion for the operator 
in both the interacting and the KS systems are written 
down, and subtracted. Since the time-dependent den- 
sity is the same in both systems, the difference vanishes. 
Usually, the Hartree term also separately satisfies the re- 
sulting equation, and so can be subtracted from both 
sides, yielding a condition on the XC potential alone. 
This procedure is well-described in the Vignale chapter 
for the zero XC force theorem. 

Zero XC force and torque: These are very simple 
conditions saying that interaction among the particles 
cannot generate a net force [30, 31]: 

J d^r n(r, t) Vwxc(r, i) = (30) 
j d\ n{r,t) r x V«xc(r,t) = J d\ r x 

where jxc(r,<) is the difference between the interacting 
current density and the KS current density [32]. The 
second condition says that there is no net XC torque, 
provided the KS and true current densities are identical. 
This is not guaranteed in TDDFT (but is in TDCDFT). 
The X-only KLI approximation, though incredibly accu- 
rate for ground state DFT, was found to violate the zero- 
force condition [33] . This is because KLI is not a solution 
to an approximate variational problem, but instead an 
approximate solution to the OEP equations. This means 
KLI also violates the virial theorem [34], which we de- 
scribe next. 

XC Powder and Virial: By applying the same 
methodology to the equation of motion for the Hamil- 
tonian, we find [35]: 

/■ 3 dn{vt) dE^c ,o1^ 

Jdr — .xc(r^) = (31) 

while another equation of motion yields the virial theo- 
rem, which intriguingly has the exact same form as in 
the ground state, Eq. (11): 

-J d\n{rt) r-Vi;xc["](rt) = E^c[n]{t)+Ta[n]{t). (32) 

These conditions are so basic that they are trivially satis- 
fied by any reasonable approximation, including ALDA, 
AA, and XX. Thus they are more useful as detailed 
checks on a propagation scheme, as mentioned earlier. 
The correlation contribution to the latter is very small. 
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and makes a very demanding test. But because the en- 
ergy does not play the same central role as in the ground- 
state problem (and the action is not simply the time- 
integral of the energy - see Robert's chapter 2), testing 
the propagation scheme is all they are used for so far. 

III. 3. Self-interaction 

For any one-electron system, 

^x(r, t)^- [ d\' Voir, t) = (N^ 1), 

J \r~r'\ 

(33) 

These conditions are automatically satisfied by XX. 
These conditions are instantaneous in time, so any adia- 
batic approximation that satisfies the ground-state con- 
ditions of Eq. (14) will also satisfy these time-dependent 
conditions, e.g. AA. On the other hand, LDA violates 
self-interaction conditions in the ground-state, so ALDA 
also violates these conditions in TDDFT. 



III. 4. Initial-state dependence 

There is a simple condition based on the principle that 
any instant along a given density history can be regarded 
as the initial moment [36, 37]. This follows very natu- 
rally from the fact that the Schrodinger equation is first 
order in time. When applied to both interacting and 
non-interacting systems, we find: 

v^c[n; ^{t'), Ht')]{rt) - v^c[n; ^-(0), $(0)](rt) for t > t' , 

(34) 

This is discussed in much detail in Nccpa's chapter. Here 
we just mention that any adiabatic approximation, by 
virtue of its lack of memory and lack of initial-state de- 
pendence, automatically satisfies it. Interestingly, al- 
though XX is instantaneous in the orbitals, it has mem- 
ory (and so initial-state dependence) as a density func- 
tional (when applied to more than two unpolarized elec- 
trons). 

This condition provides very difficult tests for any func- 
tional with memory. Consider any two evolutions of an 
interacting system, whose wavefunctions ^E* and ^E*' be- 
come equal after some time, t^. This condition requires 
that the non-interacting systems have identical XC po- 
tentials at that time and forever after, even though they 
had different histories before then. This is illustrated in 
Fig. 2. An approximate functional with memory is un- 
likely, in general, to produce such identical potentials. 

III. 5. Coupling-constant dependence 

Because of the lack of a variational principle for the en- 
ergy, there are no definite results for various limits, as in 
Eq. (13), nor is there a simple extension of the adiabatic 





^(t) 




— ^'(t) 






^ ~ — ^ 


tc t ^ 


1 


— ^^xc[^] 
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FIG. 2. An illustration of the condition based on initial state 
dependence. The two wavefunctions "if and '5' become equal 
at time tc, and therefore the KS potentials must become equal 
then and forever after. Color online. 

connection formula (7), though Gorling proposed an ana- 
log for time-dependent systems [38]. But there remains 
a simple connection between scaling and the coupling- 
constant for the XC potential [35]. For exchange, analo- 
gous to Eq. (12), the relation is linear: 

z;x[n^;$7(0)](ri) =7«x[n;$(0)](7r,72t), (35) 

where 

$^(0) = 7'^/' *(7ri , . . . , 7rAr; i - 0) (36) 

is the normalized initial state of the Kohn-Sham system 
with coordinates scaled by 7, and, for time-dependent 
densities, 

n^{r,t) = -f''n{jr,jh). (37) 

There is no simple correlation scaling, but we can relate 
the coupling-constant to scaling and find, analogous to 
Eq. (7): 

«^[n;vl/(0),$(0)](rt) = 

X'vc[n,/y, ^i/^{0), $1/^(0)] (Ar, Xh), (38) 

where ^"1/^(0) is the scaled initial state of the interacting 
system, defined as in Eq. (36) and replacing 7 with 1/A. 
For finite systems, it seems likely that taking the limit 
A — >■ makes the exchange term dominant (just as in the 
ground-state) [26], but this has yet to be proven. 

III. 6. Translational invariance 

Consider a rigid boost X(i) of a system starting in its 
ground state at t = 0, with X(0) = dX/dt{0) = 0. Then 
the exchange-correlation potential of the boosted density 
will be that of the unboosted density, evaluated at the 
boosted point, i.e., 

Vxc [n'] (r , i) = v^a [n] (r - X(t) , <) , (39) 
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Consequences of general conditions 



Linear response 



where n'(r,t) = n(r — X(i),i). This condition is univer- 
sally valid [30]. The GK approximation was found to 
violate this condition, which spurred on the development 
of the VK approximation. 



IV. LINEAR RESPONSE 

In the special case of linear response, all exchange- 
correlation information is contained in the kernel /xc- 
Linear response is utilized in the great majority of 
TDDFT calculations, and Strubbe thoroughly discusses 
the methods involved in Chapter 7. As explained in 
Chapter 24 (Martin Head-Gordon) and Ref. [24], the 
chief use of linear response has been to extract electronic 
excitations. In this section, we shall discuss the exact 
conditions that pertain to /xc, regardless of how it is 
employed. 



IV. 1. Consequences of general conditions 

Each of the conditions listed below for /xc can be de- 
rived from a general condition in Section III. 

Adiabatic limit: For any finite system, the exact 
kernel satisfies: 

lim/xc(r,r';..)= /''^-;P (40) 
w-i-o dn[r)dn[r) 

where E^q is the exact XC energy. Obviously, any adi- 
abatic functional satisfies this, with its corresponding 
ground-state approximation on the right. 

Zero force and torque: The exact conditions on the 
potential of Section III. 2 also yield conditions on /xc, 
when applied to an infinitesimal perturbation (see Vig- 
nale chapter). Taking functional derivatives of Eq. (31) 
yields 



/ 



n(r) V/xc(r,r» = -V'«xc(r') (41) 



and 



Jd\ n{r) r x V/xc(r, r'; c^) = -r' x V'i>xc(r'), (42) 

the latter assuming no XC transverse currents. Again, 
these are satisfied by ground-state DFT with the static 
XC kernel, so they are automatically satisfied by any adi- 
abatic approximation. Similarly, in the absence of cor- 
relation, they hold for XX. The general conditions cm- 
ploying energies, Eqs. (31) and (32), do not yield simple 
conditions for the kernel, because the functional deriva- 
tive of the exact time-dependent XC energy is not the 
XC potential. 



Self-interaction error: For one electron, functional 
differentiation of Eq. (33) yields: 

/x(r,r';L^) = -l/|r-r'|, /c(r,r';^)-0 (A^ - 1). 

(43) 

These conditions are trivially satisfied by XX, but vio- 
lated by the density functionals ALDA, GK, and VK. 

Initial-state dependence: The initial-state condi- 
tion, Eq. (34), leads to very interesting restrictions on 
/xc for arbitrary densities. But the information is given 
in terms of initial-state dependence, which is very diffi- 
cult to find. 

Coupling-constant dependence: The exchange 
kernel scales linearly with coordinates, as found by dif- 
ferentiating Eq. (35): 

/x[n^](r,r',cj) =7/x[n](7r,7r',cj/72). (44) 

A functional derivative and Fourier-transform of Eq. (38) 
yields [39] 



(45) 



These conditions are trivial for XX. They can be used to 
test the derivations of correlation approximations in cases 
where the coupling-constant dependence can be easily 
deduced. More often, they can be used to generate the 
coupling-constant dependence when needed, such as in 
the adiabatic connection formula of Eq. (7). 

A similar condition has also been derived for the 
coupling-constant dependence of the vector potential in 
TDCDFT [40]. 



IV. 2. Properties of the kernel 

The kernel has many additional properties that come 
from its definition and other physical considerations. 

Symmetry: Because the susceptibility is symmetric, 
so must also be the kernel: 



/xc(r,r';a;) = /xc(r',r;a;) . 



(46) 



This innocuous looking condition is satisfied by any adi- 
abatic approximation by virtue of the kernel being the 
second derivative of an energy, and is obviously satisfied 
by XX. 

Kramers— Kronig: The kernel /xc(r,r',a;) is an an- 
alytic function of to in the upper half of the complex 
w-plane and approaches a real function /xc(i", r'; oo) for 
a; — > oo. Therefore, defining the function 

A/xc(r,r',c^) = /xc(r,r',^)-/xc(r,r';oo), (47) 

we find 

~doj' 3/xc(r,r', w') 



5iA/xc(r,r',c.)=7' 



(48) 
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Properties of the kernel 



Linear response 



and 

J TT uj' — UJ 

The kernel /xc(r, r'; i, i') is real- valued in the space and 
time domain, which leads to the condition in the fre- 
quency domain, 

/xc(r,r';c.) = /*Jr,r';-c.). (50) 

The simple lesson here is that any adiabatic kernel (no 
frequency dependence) is purely real, and any kernel with 
memory has an imaginary part in the frequency domain 
(or else is not sensible). Many of the failures of cur- 
rent TDDFT approximations, e.g. the fundamental gap 
of solids, are linked to the lack of an imaginary part of the 
kernel [41]. Because adiabatic approximations produce 
real kernels, we see that memory is required to produce 
complex kernels. Hellgren et al. [42] showed that XX has 
a complex kernel, since it has frequency-dependence (for 
more than 2 electrons). Both GK and VK have complex 
kernels satisfying the Kramers-Kronig conditions. 

Adiabatic connection: A beautiful condition on the 
exact XC kernel is given simply by the adiabatic connec- 
tion formula for the ground-state correlation energy: 




/ dA5[x\r,r';^)-Xs(r,r';^)] . (51) 

Combined with the Dyson-like equation of Chapter 1 for 
as a function of Xs and /xc , this is being used to gen- 
erate new and useful approximations to the ground-state 
correlation energy [43, 44]. Although computationally 
expensive, ways are being found to speed up the calcula- 
tions [45]. 

Eq. (51) provides an obvious exact condition on any 
approximate XC kernel for any system. Thus every sys- 
tem for which the correlation energy is known can be 
used to test approximations for /xc. Note that, e.g, us- 
ing ALDA for the kernel implicit in (51) does not yield 
the corresponding E}^^^, but rather a much more sophis- 
ticated functional [39] . Even insertion of /x yields corre- 
lation contributions to all orders in Ec- And lastly, even 
the exact adiabatic approximation, /xc[?^o](r, r'; = 0), 
does not yield the exact Exc[no]- 

Functional derivatives: A TDDFT result ought to 
come from a TDDFT calculation, but this is not always 
the case. By a TDDFT calculation, we mean the result 
of an evolution of the TDKS equations of chapter 1 with 
some approximation for the XC potential that is a func- 
tional of the density. This implies that the XC kernel 
should be the functional derivative of some XC poten- 
tial, which also reduces to the ground-state potential in 
the adiabatic limit. All the approximations discussed 



here satisfy this rule. But calculations that intermix ker- 
nels with potentials in the solution of Casida's equations 
violate this condition, and run the risk of violating un- 
derlying sum-rules. 

IV. 3. Excited states 

The following conditions have to do with the challenges 
of obtaining excited states in the linear response regime. 

Infinite lifetimes of eigenstates: This may seem 
like an odd requirement. When TDDFT is applied to 
calculate a transition to an excited state, the frequency 
should be real. This is obviously true for ALDA and 
exact exchange, but not so clear when memory approx- 
imations are used. As mentioned in Section IV. 2, the 
Kramers-Kronig relations mean that memory implies 
imaginary XC kernels, and these can yield imaginary con- 
tributions to the transition frequencies. Such effects were 
seen in calculations using the VK for atomic transitions 
[46]. Indeed, very long lifetimes were found when VK 
was working well, and much shorter ones occurred when 
VK was failing badly. 

Single-pole approximation for exchange: This is 
another odd condition, in which two wrongs make some- 
thing right. Using Gbrling-Levy perturbation theory 
[47], one can calculate the exact exchange contributions 
to excited state energies [48, 49]. To recover these re- 
sults using TDDFT, one does not simply use /x, and 
solve the Dyson-like equations. Like with Eq. (51), the 
infinite iteration yields contributions to all orders in the 
coupling-constant . 

However, the single-pole approximation truncates this 
series after one iteration, and so drops all other orders. 
Thus the correct exact exchange results are recovered 
in TDDFT from the SPA solution to the linear response 
equations, and not by a full solution [50] . This procedure 
can be extended to the next order [51]. 

Double excitations and branch cuts: Maitra et 
al [52, 53] argued that a strong w-dependence in /xc al- 
lows double excitation solutions to Casida's equations, 
which effectively couples double excitations to single ex- 
citations. Similarly, the second ionization of the He atom 
implies a branch cut in its /xc at the frequency needed 
[54]. Under limited circumstances, this frequency depen- 
dence can be estimated, but a generalization [55] has been 
proposed. It would be interesting to check its compliance 
with the conditions listed in this chapter. 

Excitations in the adiabatic approximation: One 

misleading use of linear response has been to test the 
quality of different approximations to the ground-state 
i?xc- For instance, Jacquemin et al. [56] calculated 
the excitation energies for approximate E^c functionals 
within adiabatic TDDFT and compared them to exper- 
imental values. However, even within AA - using the 
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Excited states 



Extended systems and currents 



adiabatic approximation with the exact i?xc ~ the ex- 
act excitations would not be not obtained. Thus a good 
ground-state E^c used in adiabatic linear response will 
not necessarily give good excitation energies. 

Scattering theory and real-time propagation: A 

vastly under-appreciated exact condition for TDDFT is 
the equivalence of time-dependent propagation and scat- 
tering theory. This can be particularly important in un- 
derstanding the relation between bound and continuum 
states. 

For example, much early work in TDDFT was per- 
formed by Yabana and Bertsch [57], propagating ALDA 
for atoms and molecules in weak electric fields. By 
Fourier transformation of the time-dependent dipole mo- 
ment, one can extract the photoabsorption spectrum. 
The fruitfly of such calculations is benzene, with a large 
TT — >■ TT* transition at about 6.5 eV, accurately given by 
ALDA. But closer inspection shows that the LDA ion- 
ization threshold is at about 5 eV, because the LDA XC 
potential is not deep enough. Thus this transition is in 
the LDA continuum, yet its position and area are given 
reasonably well by ALDA. This is no coincidence: ALDA 
describes the time-dependent density and its propagation 
for moderate times very well. All that has changed is the 
choice of complete set of states onto which to project the 
results! 

By following this logic, Wasserman et al. [58] could 
capture the effect of Rydberg transitions using ALDA. 
However, ALDA puts many bound states in the contin- 
uum due to the exponential fall-off of the KS-LDA po- 
tential (as mentioned in Section 11). Thus the ionization 
potentials for the ALDA states are wrong, but the oscil- 
lator strength in the LDA continuum accurately approx- 
imates that of the true Rydberg transitions to the exact 
bound states. (However, it is not an exact condition that 
the KS oscillator strengths be correct, not even at the 
threshold where KS captures the right energy [59].) Us- 
ing a trick due to Fano [60], Wasserman showed [61] that 
the quantum defect, an excruciatingly sensitive measure 
of the Rydberg transition frequencies, could be extracted 
from ALDA. Ref. [62] shows the accuracy of this calcu- 
lation for He, Be, and Ne, whereas Ref. [63] shows the 
qualitative failure of ALDA for transitions to high angu- 
lar momentum eigenstates (starting at the d orbitals). 

One can go further, and even consider true continuum 
states. In scattering theory, the continuum states of the 
A^-|-l particle system describe how a single electron scat- 
ters from an N particle system. Wasserman [64] and van 
Faassen [65] developed methods to calculate scattering 
amplitudes and phase shifts based on time-propagation 
within TDDFT. With a given approximation, one can 
calculate the susceptibility of an atomic anion and de- 
duce the scattering amplitude for an incident electron 
[66]. 

Both these examples (the quantum defect and scat- 
tering) can be connected in the same framework [67], 
and they illustrate that TDDFT fundamentally concerns 



time-propagation. Present-day approximations yield 
promising results; simple approximations like ALDA of- 
ten yield accurate time-dependent densities, but their 
projection onto individual Kohn-Sham eigenstates may 
appear far more complicated. 



V. EXTENDED SYSTEMS AND CURRENTS 

As mentioned in Section 1.4, care must be taken when 
extending exact ground-state DFT results to extended 
systems. This is even more so the case for TDDFT. 
The first half of the RG theorem (chap 1) provides a 
one-to-one correspondence between potentials and cur- 
rent densities, but a surface condition must be invoked 
to produce the necessary correspondence with densities. 
Without this condition, it can readily be seen that two 
periodic systems with completely different physics can 
have the same density [68], as in Fig. 3. With hindsight, 
this is very suggestive that time-dependent functionals 
may contain a non-local dependence on the details at a 
surface. As such, they are more amenable to local ap- 
proximations in the current rather than the density. 




FIG. 3. Electrons on a ring. A magnetic field B(t) is turned 
on and steadily increases in (b); the resulting electric field 
E(r) is uniform on a thin ring, accelerating electrons around 
the ring, producing the probability current j(r, t). Note that 
in both (a) and (b) the densities are equal. Color online. 



V.l. Gradient expansion in the current 

As discussed elsewhere (Vignalc chapter) and first 
pointed out by Dobson [69], the frequency-dependent 
LDA (GK approximation) violates the translational in- 
variance condition of Section 111.6. One can trace this 
failure back to the non-locality of the XC functional in 
TDDFT. But, by going to a current formulation, every- 
thing once again becomes reasonable. The gradient ex- 
pansion in the current, for a slowly varying gas, was first 
derived by Vignale and Kohn [70], and later simplified 
by Vignale, Ullrich, and Conti [71], and is discussed in 
much detail in the Vignale chapter. 
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For our purposes, the most important point is that, by 
construction, VK satisfies translational invariance. The 
frequency-dependence shuts off' (it reduces to ALDA) 
when the motion is a rigid translation, but turns on when 
there is a true (non-translational) motion of the density 
[70]. 

Any functional with memory should recover the VK 
gradient expansion in this limit, or justify why it does 
not. However, the VK approximation is only the gra- 
dient expansion, which for the ground-state was found 
to violate sum rules, as mentioned in Section 1.4. It is 
therefore likely that there exists something like a gen- 
eralized gradient approximation, which is more accurate 
than VK. 



V.2. Polarization of solids 

A decade ago, GGG [72] pointed out that the periodic 
density in an insulating solid in an electric field is insuf- 
ficient to determine the one-body potential, in apparent 
violation of the Hohenberg-Kohn theorem [25]. In fact, 
this effect appears straightforwardly in the static limit 
of TDCDFT, and is even estimated by calculations using 
the VK approximation [68, 73]. When translated back to 
TDDFT language, one finds a l/g^ dependence in /xc, 
where q is the wavevector corresponding to r — r'. This 
requires /xc to have the same degree of nonlocality as 
the Hartree kernel, and this is missed by any local or 
semilocal approximation, such as ALDA, but is built in 
to XX [74] or AA. The need for a. contribution in 
the optical response of solids led to much development 
[75] for a kernel that allows excitons [76, 77]. Since the 
RG theorem can be proven for solids in electric fields of 
nonzero q, one can extract the g — >■ (a constant E field) 
result at the end of the calculation [68] . 

VI. SUMMARY 

What lessons can we take away from this brief survey? 



1. In the ground-state theory, the total XC energy is 
crucial for determining the energy of the system, 
and many conditions are proven for that functional. 
This is not so for TDDFT, for which only the time- 
dependent density matters. In the non-interacting 
system, the KS potential, and specifically its XC 
component, is what counts. 

2. Explicit density functionals have poor-quality po- 
tentials, e.g. LDA and GGA. Thus successes in 
ground-state DFT do not translate directly into 
successes in TDDFT. One of the greatest challenges 
is that the potential is a far more sensitive func- 
tional of the density than vice versa. Though we 
have enumerated many conditions on the XC po- 
tential, it is important to determine which condi- 
tions significantly affect the density, including those 
aspects of the density that are relevant to experi- 
mental measurements. 

3. The adiabatic approximation satisfies many exact 
conditions by virtue of its lack of memory. Inclu- 
sion of memory may lead to violations of condi- 
tions that adiabatic approximations satisfy. This is 
reminiscent of the ground-state problem, where the 
gradient expansion approximation violates several 
key sum rules respected by the local approximation. 
Explicit imposition of those rules led to the devel- 
opment of generalized gradient approximations. 

As shown in several chapters in this book, many people 
are presently testing the limits of our simple approxima- 
tions, and very likely, these or other exact conditions will 
provide guidance on how to go beyond them. 
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